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ABSTRACT 

The  interaction  between  the  tire  and  terrain  has  long  been  of  interest  for  vehicle  dynamic  simulation.  Detailed  tire 
models  produce  accurate  results,  but  are  too  computationally  intensive  for  the  iterative  vehicle  design  process.  The 
objective  of  this  work  is  to  develop  a  computationally  efficient  way  to  estimate  the  tire’s  deformed  shape,  from 
which  the  tire  forces  can  be  evaluated.  A  novel,  planar,  quasi-static,  constraint  mode  tire  model  is  developed  to 
address  this  objective.  The  required  model  parameters  are  reduced  to  two  non-dimensional  shape  parameters  and  an 
overall  stiffness  factor.  An  admissible  domain  of  the  shape  parameters  is  developed  based  on  the  deformation 
limitations  of  a  physical  tire.  Specifically,  no  single  harmonic  may  dominates  the  tire  shape  and  the  low  spatial 
frequency  components  must  contribute  more  than  higher  frequency  components  to  the  overall  tire  shape.  The  ability 
of  the  model  to  accurately  predict  the  spindle  force  is  evaluated  by  comparing  simulation  and  experimental 
responses  for  quasi-static  flat  plate  and  cleat  tests.  This  work  provides  a  simple,  accurate  tire  model  for 
circumferential  displacement  and  vertical  spindle  force  prediction  to  improve  the  vehicle  design  process. 
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NOMENCLATURE 


N  Total  number  of  tire  segments 
u  Radial  deflection. 

R  Tire  radius. 

F„  Radial  force  acting  on  tire  segment. 

F  Vertical  spindle  force. 
f  Linearly  distributed  load  density. 

K  Radial  stiffness  density. 

Ui,  Bending  potential  energy  for  each  tire  segment. 
M  Bending  moment  for  each  tire  segment. 

E  Elastic  modulus. 

/  Second  moment  of  inertia  for  each  tire  segment, 
l/j  Shear  potential  energy  for  each  tire  segment. 

G  Shear  modulus. 

A  Cross-sectional  area. 

Ug  Elastic  potential  energy  for  each  tire  segment. 
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Distributed  physical  bending  stiffness, 

72  Distributed  physical  shear  stiffness. 

73  Distributed  physical  radial  stiffness. 
kg  Overall  model  stiffness. 

a^  a2  Model  shape  parameters. 


INTRODUCTION 

With  the  advent  of  autonomous  vehicles  and  the  simultaneous  pressure  to  improve  energy  efficiency,  there  are  ever 
increasing  demands  on  vehicle  engineers.  Clearly  there  are  more  electronic  and  mechatronic  systems,  more  sensors 
and  more  actuators.  Batteries  increase  vehicle  weight  at  the  same  time  that  the  incentive  to  improve  fuel  economy 
continues  to  grow,  which  in  turn  drives  a  reduction  in  the  overall  vehicle  weight.  These  new  technologies  are 
making  vast  improvements  in  vehicle  performance  and  the  way  drivers  think  of  their  vehicles.  They  also  introduce 
new  challenges  to  vehicle  engineers. 

The  chassis  engineer  is  faced  with  numerous  design  challenges  (weight,  cost,  ride,  handling,  noise,  reliability, 
packaging...)  while  interacting  with  several  influential  groups  (scientific  labs,  proving  ground,  suppliers,  safely, 
manufacturing,  powertrain  and  body  groups)  that  result  in  a  highly  iterative  design  process.  All  these  challenges  and 
groups  must  be  satisfied,  to  the  best  of  the  engineer’s  ability,  within  a  very  short  timeframe.  Within  this  cacophony 
of  ever-changing  requirements  and  influential  groups,  some  fundamental  physical  truths  remain.  Regardless  of  how 
the  power  is  developed,  or  how  the  wheels  are  steered,  all  the  vehicle  forces  must  ultimately  be  reacted  through  the 
vehicle  suspension  and  tires.  The  focus  of  this  work  is  the  interaction  between  the  tires  and  the  main  excitation  to 
the  vehicle:  the  terrain  [1]. 

This  work  addresses  the  means  by  which  non-deformable  terrain  imposes  a  unilateral  geometric  boundary  constraint 
on  rolling  tires  to  which  the  chassis  responds  by  generating  loads,  moments,  motions,  deformations,  etc.  The  terrain 
surface  remains  a  consistent  excitation  to  the  vehicle  throughout  the  iterative  design  process.  Understanding  and 
properly  modeling  this  interaction  is  critical  in  predicting  the  loads  and  the  resulting  vehicle  performance.  However, 
actual  loads  are  only  available  at  the  conclusion  of  the  design  and  development  process  when  it  is  difficult  and 
expensive  to  make  modifications. 

Clearly,  computationally  efficient  and  accurate  vehicle  dynamic  simulations  are  critical  throughout  this  iterative 
process.  High  fidelity,  but  computationally  intensive  models  that  provide  very  accurate  information  after  the  parts 
have  been  released  are  of  no  use.  Simulation  results  arising  from  overly  simplified  models  that  cannot  be  trusted 
and  have  little  value.  This  work  aims  to  strike  a  balance  between  heuristic  tire  models  (such  as  a  linear  point- 
follower)  that  lack  the  fidelity  to  make  accurate  chassis  load  predictions  and  computationally  intensive  models  that 
cannot  provide  timely  predictions. 

The  objective  of  this  work  is  to  develop  a  computationally  efficient,  planar  tire  model  that  accurately  predicts  the 
lower-frequency,  but  not  necessarily  low  amplitude,  tire  shape.  The  emphasis  is  intentionally  placed  on  the  shape  of 
the  deformed  tire,  which  must  be  capable  of  representing  the  bridging  property  (over  narrow  cracks  in  the  pavement) 
and  enveloping  property  (over  sharp  bumps).  This  scope  includes  relatively  large  deformations,  approaching,  but 
not  including,  the  point  of  rim  strike  (when  the  tire  sidewall  is  folded  such  that  the  upper  portion  is  in  direct  contact 
with  the  lower  portion  and  there  is  a  discontinuous  and  dramatic  increase  in  apparent  stiffness).  A  simple  linear 
relationship  between  the  tire  circumferential  deformation  and  the  spindle  load  prediction  is  used  to  develop  the 
deformed  tire  shape  with  the  understanding  that  more  complex  force-deflection  relationships  could  be  developed 
using  this  deformed  shape  as  an  initial  approximation. 

Higher-frequency  dynamics  such  as  those  that  would  be  required  for  noise  predictions  are  outside  the  scope  of  this 
work.  The  focus  is  the  deformed  shape  of  the  tire  in  the  region  of  the  contact  patch.  The  majority  of  the  tire  is  non¬ 
contacting  and  must  adopt  a  smooth  shape.  That  is,  the  tire  model  parameterization  must  be  such  that  the  deformed 
tire  shape  does  not  exhibit  higher  frequency  dynamics  that  are  typically  associated  with  higher  free-vibration  normal 
modes.  Defining  the  parameter  constraints  such  that  these  higher  order  dynamics  are  omitted  is  one  of  the  major 
contributions  of  this  work.  This  results  in  the  development  of  an  admissible  region  for  the  tire  shape  parameters. 
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Developing  a  parsimonious  parameter  set  is  one  of  the  major  contributions  of  this  work.  The  tire  model  developed 
in  this  work  is  completely  parameterized  by  three  parameters:  a  single  stiffness  parameter  and  two  shape  parameters. 
Experimental  data  is  used  to  validate  the  approach  using  two  quasi-static  tests,  a  flat  plate  and  a  cleat  test,  with 
relatively  large  deformations.  The  model  is  capable  of  bridging  or  enveloping  small  irregularities  in  the  terrain 
surface  and  generating  the  correct  tire  shape  over  irregular  terrain. 

The  remainder  of  this  work  is  developed  as  follows.  Some  background  on  planar  tire  modeling  and  the  use  of 
constraint  modes  is  briefly  reviewed.  The  mathematical  derivation  of  the  constraint  mode  model  is  developed  such 
that  the  entire  model  can  be  described  by  just  three  parameters.  The  constraints  on  the  tire  shape  parameters  are 
developed  based  on  the  required  physical  properties  of  a  tire  and  the  requirement  that  the  shape  of  the  non¬ 
contacting  region  of  the  tire  must  be  smooth.  Finally,  the  parameterization  process  for  a  specific  tire  is  developed 
and  validated  with  quasi-static  experimental  results. 

BACKGROUND 

A  considerable  amount  of  tire  modeling  research  has  been  conducted  and  a  comprehensive  review  of  the  current 
state  of  the  art  is  presented  by  Umsrithong  [2].  Although  a  complete  review  of  recent  developments  is  not  presented 
herein,  the  results  from  several  of  these  works  have  direct  relevance  to  the  objectives  of  this  study.  For  example, 
Badalamenti  et  al.  showed  that  a  radial  spring  tire  model,  in  which  the  radial  spring  element  deflection  depends  on 
the  adjacent  element  deflections,  could  accurately  describe  the  tire  enveloping  behavior  in  an  efficient  model  [3]. 
Presently,  a  planar  ring  model  with  an  elastic  foundation  that  emphasizes  the  circumferential  deformation  is 
developed  to  accurately  and  efficiently  describe  the  tire  enveloping  behavior. 

The  deformed  shape  of  the  tire  (including  bridging  and  enveloping)  was  accurately  predicted  by  Zegelaar  and 
Pacejka  by  using  a  flexible  ring  to  simulate  the  quasi-static  response  of  a  tire  rolling  over  an  uneven  surface  [4]. 
Although  the  simulation  results  agreed  with  the  experimental  data,  the  complexity  of  the  model  requires  additional 
computational  requirements.  The  objective  to  verify  the  parameter  identification  process  via  the  deformed  shape  of 
the  tire  is  accomplished  using  experimental,  single  point-load  testing.  This  single  point  excitation  to  the  tire  model 
is  also  used  analyze  the  frequency  content  of  the  deformed  shape  of  the  tire.  Presently,  single  point-load 
experimental  results  are  used  to  parameterize  the  constraint  mode  tire  model  developed  in  this  work.  A  similar 
approach  was  developed  by  Loo,  who  modeled  the  tire  by  using  a  more  complex  flexible  circular  ring  under  tension 
with  a  nest  of  linear  springs  and  dampers  arranged  radially  and  experimentally  identified  the  parameters  using  single 
point-load  tests  [5].  However  in  Loo’s  study  the  ring  tension  and  foundation  stiffness  requires  contact  patch 
measurements. 

It  is  assumed  in  this  work  that  the  low-frequency  deformation  of  the  tire  shape  can  be  adequately  defined  by  the 
quasi-static  constraint  modes.  This  assumption  is  given  some  credence  by  Gillespie’s  study  in  which  a  radial  spring 
model  is  used  to  simulate  the  stiffness  variation  circumferentially  and  the  magnitude  of  the  radial  force  variation  is 
found  to  be  relatively  independent  of  speed  [6].  It  is  further  assumed  in  this  work  that  the  transient  response  can  be 
adequately  represented  with  a  mass-spring  tire  model.  Takayama  et  al.  develops  such  a  model  to  predict  the 
transient  response  of  a  tire  encountering  a  cleat,  where  the  belt  and  tread  region  is  modeled  by  a  rigid  ring  and 
deflections  from  the  cleat  are  absorbed  by  a  linear  and  planar  spring  attached  to  the  rigid  ring  [7]. 

Ferris  demonstrates  that  a  static  constraint  mode  tire  model  can  be  used  to  capture  the  tire  enveloping  and  bridging 
properties  with  an  axisymmetric  and  circumferentially  isotropic  model  [8].  It  is  modeled  as  an  inextensible,  but 
flexible,  ring  supported  by  an  elastic  foundation  represented  by  the  radial  stiffness,  k,  as  shown  in  Figure  la.  The 
tire  is  divided  into  A  segments  and  each  segment  is  considered  an  Euler  elastic  beam  of  length  ItiR/N  and  cross- 
sectional  area.  A,  shown  as  a  shaded  region  in  Figure  lb.  The  radial  deflection  vector  {«}  is  used  to  describe  the 
circumferential  tire  displacement  of  each  tire  segment.  The  quasi-static  constraint  modes  are  developed  from 
Equation  (1). 
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Figure  1.  Schematic  plot  of  constraint  tire  model  and  a  piece  of  tire  segment 


[A-llu]  =  {F}  (1) 

The  physical  degrees  of  freedom  are  categorized  into  active  and  omitted  sets,  represented  by  superscripts  ‘a’  and  ‘o’ 
respectively.  In  general,  the  boundary  degrees  of  freedom  must  be  a  subset  of  the  active  degrees  of  freedom.  In  the 
tire  modeling  case,  the  active  degrees  of  freedom  must  include  those  in  contact  with  the  terrain.  Therefore  the  mass 
and  stiffness  matrices  can  be  reordered  and  partitioned  as  given  in  equation  (2). 

.[^aa]  r(F“}| 

[«■““]  [«■““]]  I  {0)1 

An  iterative  process  to  identify  the  active  constraints  at  the  tire-terrain  interface  is  formulated,  as  shown  by  Figure  2, 
beginning  with  the  undeformed  tire  shape  shown  in  Figure  2(a).  The  segment  with  the  largest  interference  is 
considered  the  first  active  constraint,  {u“}  shown  as  a  single  small  circle  in  Figure  2(b).  The  omitted  constraint 
vector  (u")  can  be  obtained  as  Equation  3  through  Guyan  mode  reduction  [9]. 

{u“}  =  (3) 

The  resulting  circumferential  deformation  with  the  first  active  constraint  is  predicted,  as  shown  in  Figure  2(b)  where 
the  active  constraint  is  indicated  by  a  small  circle.  The  segment  with  the  most  interference  in  Figure  2(b)  is  selected 
as  the  second  active  constraint  and  a  new  tire  shape  is  computed  as  shown  in  Figure  2(c).  This  process  proceeds 
until  there  is  negligible  interference,  as  shown  by  Figure  2(f). 


(a) 


(c) 


(d)  (e)  (I) 

Figure  2.  The  iterative  process  to  identify  the  active  constraints 

The  objective  of  this  work  is  to  properly  define  this  stiffness  matrix  [F].  The  stiffness  matrix  is  defined  with  a  few 
parameters  (i.e.,  one  stiffness  and  two  shape  parameters)  and  gives  reasonable  shape  predictions  for  relatively  large 
deflections.  Perhaps  most  importantly,  the  resulting  model  yields  shapes  that  are  intuitive;  the  tire  shape  follows  the 
contour  of  the  terrain  along  the  contact  area  and  is  smooth  in  the  non-contacting  region  of  the  tire. 

MODEL  DEVELOPMENT 

The  quasi-static  equation  resulting  from  the  application  of  Hamilton’s  principle  for  one  tire  segment  is  written  as 
Equation  (4). 

rt2  (4) 

(5U^  -  SU^  -  5U,  +  SU^)dt  =  0  Vti, 

■'ti 
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The  virtual  changes  are  computed  as  Equations  (4)  -  (8). 


2nR  d^u 

SUu  = - El  ——Su 

"  N  dx^ 

2nR 

SUp  — - KuSu 

N 

2nR  d^u 
SUs  ^  — —  GA  Su 

N  dx^ 

2nR 

SVyp  = 


The  equation  of  motion  for  each  segment  is  then  simplified  to  Equation  (9). 

d^u 


(5) 

(6) 

(7) 

(8) 

(9) 


For  simplicity,  the  following  parameters  are  defined  in  Table  1.  The  tire  stiffness  parameters  (yi,y2and  Ys)  ^re 
invariant  for  a  specific  tire;  the  stiffness  associated  with  a  particular  tire  segment  is  a  function  of  the  number  of  tire 
segments,  N.  It  should  be  clear  that  each  of  the  three  physical  stiffness  parameters  must  be  positive. 


Table  1 .  Physical  stiffness  parameters  as  a  function  of  physical  properties 


Yr  =  EI{2nRy^ 

Segment  bending  stiffness 

Y2  =  GA(2nRy^ 

Segment  shear  stiffness 

Y:3  =  K(2nR) 

Segment  radial  stiffness 

F„=^(27rR) 

Radial  force  on  n‘^  tire  segment 

The  fourth  spatial  derivative  of  u  for  the  tire  segment  is  approximated  by  adjacent  points  via  a  finite  difference 
method  as  given  in  Equation  (10). 


Uji—2  4Uj2— 1  T  6Uj2  4Uj2+i  T  trj2+2  (id) 

2nR^* 

N  ) 

The  equation  of  motion  for  the  tire  segment  is  expressed  as  Equation  (11), 


kgCt2Un-2  "b  +  kgU-n  +  kgCtiU-n^i  +  fco®2^n+2  ~ 


(11) 


where  the  stiffness  and  shape  parameters  are  defined  in  terms  of  the  tire  parameters  as  Equation  (12)  [10]. 

fco  =  6y^N^  -  2Y2N  +  Y3N-\ 

fcoCti  =  Y2N  -  4ylAf^  (12) 

fco«2  =  YiN^ 


The  symmetry  of  a  physical  tire  results  in  a  circulant  stiffness  matrix  [/f]  that  comprises  two  parts:  kg  defines  the 
overall  stiffness  and  shape  matrix  [a]  defines  the  relative  deformations  (i.e.,  the  shape)  of  the  deformed  tire  [10],  as 


shown  in  Equation  (13). 

-  1 

a-i 

0.2 

0  • 

0.2 

a-Y 

(13) 

a-i 

1 

(X2 

••  0 

(X2 

[K] 

—  kg 

(X2 

0 

a-i 

(X2 

1 

1  • 

i  0 

0 

0 

0 

=  ko{a\  =  koCirc({ai}) 

(X2 

0 

0  • 

•• 

1 . 
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In  general,  a  circulant  matrix  is  denoted  by  circ  and  eharacterized  by  the  first  row  of  the  matrix,  written  as  the 
eharacteristic  row  veetor  (aj).  The  shape  matrix  [a]  has  a  circulant  form  that  is  only  a  function  of  two  shape 
parameters  (ctj^  and  a2)  so  that  the  characterization  vector  becomes 

(aJ=[“o  %  ...  =  [1  a-^  ^2  0  ...  0]  (14) 

ADMISSIBLE  REGION  FOR  SHAPE  PARAMETERS 

Clearly  the  physical  tire  stiffnesses  must  be  positive  and  the  constraints  on  the  physical  stiffnesses  must  be  enforced 
on  the  model  parameters.  It  should  also  be  intuitive  that  the  non-contacting  region  of  a  tire  under  quasi-static 
conditions  must  adopt  a  smooth,  unrippled  shape.  That  is,  the  upper  section  of  the  tire  remains  nearly  round.  More 
specifically,  changes  in  the  radial  deformation  of  the  tire  occur  smoothly  and  in  small  increments  in  the 
circumferential  direction.  This  implies  that  the  contributions  of  higher  spatial  frequencies  must  remain  small.  This 
is  in  contrast  to  the  deformed  shape  that  a  tire  assumes  when  it  is  being  harmonically  excited,  particularly  in  the  free 
vibration  mode.  The  remainder  of  this  section  proceeds  as  follows.  First  the  physical  restrictions  on  the  tire’s 
stiffness  parameters  are  used  to  develop  a  set  of  constraints  on  the  shape  parameters.  Next,  a  single  point  load  is 
considered  in  order  to  develop  the  spatial  frequency  constraints.  The  low-frequency  spatial  constraint  is  formulated 
and  the  resulting  constraints  on  the  shape  parameters  are  developed.  This  set  of  constraints  forms  the  admissible 
region  for  the  shape  parameters. 


Physical  Constraints 


The  tire  shape  is  constrained  by  its  physical  properties  which  must  be  enforced  as  constraints  on  the  shape 
parameters  a-i  and  a2.  An  admissible  region  of  the  shape  parameter  set  is  derived  from  these  constraints. 
Specifically,  the  tire  must  possess  positive  bending,  shear  and  radial  stiffness  so  that  the  tire  stiffness  parameters, 
Y3^~^  (as  shown  in  Table  1),  must  be  positive.  The  positive  sense  for  the  deflections  and 
corresponding  forces  are  defined  such  that  the  overall  tire  stiffness,  ko  must  also  be  positive.  Therefore,  the 
relationships  developed  in  Equation  (12)  are  rearranged  to  develop  constraints  on  the  tire  shape  parameters,  a-^  and 
a2,  as  given  in  Equation  (15). 


0.2 


YiN^ 

kg 


>  0 


Y2N 

a,  +  4an  = - >  0 

ko 


(15) 


Y3N  ^ 

2ai  +  202  -I- 1  =  — ; - >  0 

ko 


Next  an  additional  constraint  that  enforces  the  low  spatial  frequency  requirement  (smoothness  for  the  non-contacting 
region  of  the  tire)  is  formulated  by  considering  the  effects  of  a  single  point  load. 


Quasi-Static  Point  Loading 


Consider  a  single  point  load  on  the  tire  such  that  only  one  degree  of  freedom  is  active.  Without  loss  of  generality, 
the  indices  of  the  degrees  of  freedom  are  assigned  such  that  this  active  point  is  the  first  degree  of  freedom.  The 
deformed  tire  shape  is  simply  computed  as  Equation  (16). 


w=(“:}=™--r;}=w 


0 


(16) 


Since  there  is  only  one  active  segment,  u“  and  feg  are  scalar  values  and  the  deformed  shape  of  the  tire,  (u)  = 
gj,  is  the  product  of  the  first  column  of  the  matrix  and  the  scalar  .  To  compute  the  inverse  matrix, 

[«]■^  an  eigen-decomposition  method  is  applied  to  the  original  shape  matrix  [a].  The  eigenvalues  of  the  shape 
matrix  [a]  are  written  in  vector  form  as  [A]  and  are  calculated  by  applying  the  Discrete  Fourier  Transform  to  the 
vector  {a}  as  shown  in  Equation  (17)  [11]. 

u}  =  Tm) 
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iV-1 


4  =  y 


k  =  0,1,2,  1 


2n  . 

W  =  e~~N^ 

The  eigenvector  corresponding  to  the  k"'  eigenvalue,  is  given  in  Equation  (18). 

1 

Vfc=-— I4^2fc  ...  ]^(N-i)fc)'  fc  =  0,1,2, 

yjN 


1 


(17) 


(18) 


The  eigenvector  matrix,  [l^],  is  defined  as  the  sequential  concatenation  of  eigenvectors,  provided  in  Equation  (19) . 

[l^]  =  [Vo  Vi  V2  •••  (19) 

The  circulant  matrix  [a]  can  be  diagonalized  by  [l^]  as  shown  in  Equation  (20),  where  [l^]^  is  the  Hermitian  matrix 

ofm. 

[a]^  [Vf  diag(Ao  -  4-i)[l^]  (20) 

The  inverse  of  [a]  is  developed  as  shown  in  Equation  (21). 

[ar^  =  [V]"diag(^  -  t^)[V] 

'to  'tjv-i 

The  inverse  matrix  of  [a]  is  also  a  circulant  matrix  because  it  shares  the  eigenvectors  of  [a]  and  the  eigenvalue 
vector  of  [a]  can  be  written  as  .  Since  the  characterization  vector  {a}  is  sparse,  it  is  convenient  to  write  the 

eigenvalues  of  [a]  as  shown  in  Equation  (22),  where  it  can  be  seen  that  the  eigenvalue  vector  {T}  is  symmetric  (i.e., 
—  ^N-k’  for  k  =  1,2, ... ,  Af  —  1). 


/2nk\  /47r/c\  (22) 

Afc  =  1  +  2aiCos  ( j  +  2a2Cos  ( j  ^  =  0,1,2  ...  Af  —  1 

The  characteristic  vector  of  the  inverse  matrix,  i.e.  the  first  row  of  is  denoted  as  .  It  can  be  obtained  by 

performing  an  Inverse  Discrete  Fourier  Transform  to  the  eigenvalue  vector  of  .  The  inverse  matrix  also 
has  a  circulant  form  as  shown  by  Equation  (23),  and  is  written  as  B  for  convenience. 


[a]  1  =  B  =  Circ({^J) 
iP]  =  W0.P1.P2 . Pn-2.Pn-i]=0^-\W7"}) 


N-1 


(23) 


The  characteristic  vector  W]  is  symmetric  (i.e.,  4  =  Pn-k.f<^  —  1-2, ...,  Af  —  1)  due  to  the  symmetry  of  {/!}•  Given 
Equation  (16)  and  the  symmetric  property  of  [p],  the  resultant  deformed  tire  shape  with  one  active  segment  is 
computed  as  Equation  (24). 
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To  summarize,  the  deformed  tire  shape  vector  [u]  for  a  single  active  element  is  the  vector  W)  scaled  by 

The  characteristic  vector  is  obtained  by  calculating  the  Inverse  Discrete  Fourier  Transformation  which 

is  determined  by  the  shape  coefficients  and  a2,  as  shown  in  Equation  (17). 
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Spatial  Frequency  Constraints 

There  are  two  spatial  frequency  constraints  placed  on  the  tire  shape.  The  first  is  that  there  is  no  single  harmonic  that 
dominates  the  shape.  The  second  is  that  low  order  (low  spatial  frequency)  components  must  contribute  more  than 
higher  order  components  to  the  overall  tire  shape.  Each  of  these  constraints  is  developed  in  turn.  First,  it  is  clear 
from  Equation  (23)  that  if  one  of  the  eigenvalues  of  the  shape  matrix,  say  /l„,  approaches  zero,  then  the 
characteristic  vector  would  be  dominated  by  the  nearly  infinite  term  and  the  remaining  components  would 
be  negligible.  In  this  case,  the  deformed  tire  shape  vector  {u}  would  be  harmonic,  resulting  in  an  unrealistic 
deformed  tire  shape.  Figure  3  shows  examples  of  undesirable  tire  shapes  that  are  dominated  by  a  single  harmonic 
term;  the  solid  line  represents  the  undeformed  tire  shape  and  the  dashed  lines  represent  the  deformed  tire  shape  with 
a  dominant  harmonic  term.  The  first  diagram,  m  =  Q,  shows  the  shape  of  a  tire  containing  only  a  constant  radial 
deformation  (which  violates  the  condition  that  perimeter  is  fixed  -  it  is  a  steel  band).  The  second  diagram  shows  the 
condition  for  m  =  3  in  which  the  bottom  of  the  tire  is  fiat,  but  the  upper-right  and  upper-left  sections  of  the  tire  are 
also  deflected  inward,  which  is  unreasonable.  These  shapes  may  correspond  to  components  of  free-vibration  of  a 
tire,  but  they  fail  to  capture  the  deformed  shape  of  a  tire  in  quasi-static  contact  with  the  pavement.  Therefore,  since 
no  single  harmonic  can  dominat  the  shape,  none  of  the  eigenvalues  of  the  shape  matrix  can  approach  zero;  this  first 
spatial  frequency  constraint  is  written  as  Equation  (25) . 

Vt  £{0,1,2  ...n-l]  (25) 


Figure  3.  Solid  line:  Undeformed  tire  shape;  Dashed  line:  deformed  tire  shape  with  dominant  harmonic 
It  can  be  shown  from  Equation  (22)  that  Tj  can  be  rewritten  as  a  ftmction  of  a  new  variable  f  as  Equation  (26). 

/lj(0  =4a2f2-l-2aif-l-l-2a2  i  =  0,1,2  ...Af  -  1  (26) 


COS 


(2nk\ 


-1,1] 


Therefore,  the  constraint  given  in  Equation  (25)  can  be  rewritten  in  terms  of  f  as  Equation  (27). 

T^(f)  ^0,Vf  £  (-1,1) 


(27) 


Recall  that  the  model  is  quasi-static  and  therefore  the  radial  deformation  vector  should  be  dominated  by  a 
combination  of  low-order  terms.  Again,  the  deformation  should  not  be  dominated  by  a  single  harmonic.  This  low- 
order  dominance  is  developed  in  mathematical  terms  by  considering  Equation  (24)  and  noting  that  low-order  tire 
shapes  occur  when  the  vector  {(]]  is  dominated  by  low-order  terms.  Next  consider  Equation  (23)  and  note  that  the 
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elements  pk  that  compose  the  vector  [p]  are  determined  by  the  relative  weighting  term,  of  the  function  W  = 

2nik .  ^ 

e  N  J ,  Therefore,  the  weighting  term  —  should  be  relatively  large  (AjCf)  should  be  small)  in  the  low  frequency 

range  (i.e.,  k  «  N/2).  This  low  frequency  range  corresponds  to  f  =  1  as  given  by  Equation  (26).  The  weighting 
term  should  be  relatively  small  (/IjCf)  should  be  large)  in  the  high  frequency  range  (i.e.,  k  »  N/2)  where  f  =  — 1. 
This  second  spatial  frequency  condition,  the  low  frequency  domination  constraint,  is  therefore  expressed  as  the 
inequality  given  in  Equation  (28)Error!  Reference  source  not  found.. 

T,(1)<T,(-1)  (28) 


Finally,  note  that  the  function  /li(f)  is  a  parabola.  Applying  the  constraint  on  a2  that  was  derived  in  Equation  (15), 
the  constraint  on  the  second  derivative  of/lj(f)  is  given  in  Equation  (29). 


=  8a2 


>  0 


(29) 


That  is,  the  parabolic  function  Tj  (f)  faces  upwards.  If  the  discriminant  is  negative,  i.e.  d  =  Aaf  —  160:2(1  — 
20:2)  <  0,  the  parabolic  function  Tj  (f)  does  not  have  any  roots.  If  the  discriminant  is  positive  there  are  two  roots 
and  if  it  is  equal  to  zero  then  there  is  one  root.  The  three  cases  are  developed  as  follows. 


Case  1.  Negative  Discriminant  A  =  4al  —  160:2(1  —  2a2)  <  0 


In  this  case,  the  function Tj (f)  has  no  roots  and  /li(f)  >  0,  Vf  £  [—1,1].  Therefore  the  constraint  shown  in 
Equation  (25)  must  hold  for  all  possible  choices  of  k  and  N.  Pairs  of  values  of  a  /  and  0:2  that  yield  a  negative 
determinant  are  sufficient  to  ensure  that  there  is  no  single  dominant  harmonic  in  the  quasi-static  response. 
Specifically,  the  inequality  in  Equation  (30)  is  a  sufficient  constraint. 


+ 


16 


<  1 


(30) 


Since  /li(0  >  0,  Vf  £  [—1,1]  and  given  the  constraint  that  a2  >  0  (from  Equation  (15)),  the  low  frequency 
domination  constraint  is  derived  from  Equation  (26)  as  Equation  (31). 

%  <  0  (31) 


Case  2.  Positive  Discriminant  A  =  Aaf  —  160:2(1  —  20:2)  >  0 

In  this  case  the  function  /li(f)  has  two  roots  given  by  Equation  (32)  where  ^2  >  fi  since  <  0  and  0:2  >  0. 

cti  -I-  Va 

“  40:2  (32) 

cti  —  Va 

The  function  /li(0  for  this  case  is  shown  in  Figure  4  as  an  upward  facing  parabola.  If  no  single  harmonic  is  to 
dominate  the  tire  shape,  then  according  to  Equation  (27)  there  must  be  no  roots  within  the  domain  of  f .  Therefore, 
the  entire  domain  off  £  (— 1,1)  must  fall  into  one  of  three  regions,  as  shown  in  Figure  4. 
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RegiOT  1 


Region  2 


Region  3 


*  \  I  '  I  /  *■ 


Figure  4.  The  function  Tj(f)  with  discriminant  A  >  0 


In  the  first  region,  the  entire  domain  of  f  £  (—1,1)  must  be  less  than  so  that  the  inequality  given  in  Equation  (33) 
must  hold. 


1  <  fi  <  ^2 

It  follows  from  Equation  (32)  that  the  inequality  expressed  as  Equation  (34)  must  also  be  true. 

+  4a2  <  — Va<  0 


(33) 


(34) 


However,  the  inequality  in  Equation  (34)  contradicts  the  physical  constraints  expressed  by  Equation  (15).  Therefore 
Region  1  is  not  admissible. 

In  the  second  region,  the  entire  range  of  Tj(f)  must  be  negative  so  that  the  inequality  given  in  Equation  (35) 
must  always  hold. 

4(-l)  <  0  (35) 


It  follows  from  Equation  (26)  that  this  inequality  is  equivalent  to  Equation  (36). 

2ai  +  202  +  1  <  0  (36) 


The  inequality  in  (36)  contradicts  the  physical  constraint  expressed  in  (15).  Therefore  Region  2  is  not  admissible. 

Finally,  in  the  third  region,  the  entire  range  of  2,j(f)  is  monotonically  increasing,  so  that  the  inequality  given  in 
Equation  (33)  must  hold. 

0<4(-1)<4(1)  (37) 


This  inequality  contradicts  the  low  frequency  domination  constraint  and  therefore  Region  3  is  not  valid.  In 
summary,  none  of  the  three  regions  for  the  positive  discriminant  are  admissible. 

Case  3.  Zero  Discriminant  A  =  —  1602(1  —  202)  =  0 

Similar  arguments  can  be  made  for  the  zero  discriminant  case  as  those  made  for  the  positive  discriminant  case.  Here 
there  is  one  root  so  that  and  ^re  effectively  equal  and  there  are  effectively  two  regions  in  Figure  4  to  consider: 
Region  1  and  Region  3.  The  conclusions  for  both  of  the  regions  are  the  same  as  in  the  case  when  the  discriminant  is 
positive.  Therefore,  the  only  case  in  which  sufficient  conditions  to  satisfy  both  spatial  frequency  constraints  is  when 
the  discriminant  is  negative  wherein  the  function  /lj(f)  has  no  roots  and  /lj(f)  >  0,  Vf  G  [—1,1]. 
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Summary  of  Admissibility  Constraints 

All  the  active  constraints  on  the  shape  coefficients  and  a2  obtained  under  the  three  physical  constraints  and  the 
two  conditions  derived  from  the  spatial  frequency  constraints  are  expressed  in  Equation  (38). 

Maf  -  16a2(l  -  20:2)  <  0  (38) 

1  <  0 

j  +  4a2  >  0 


Note  that  two  of  the  five  constraints  are  redundant.  The  resultant  admissible  region  for  and  a2 's  plotted  as 
Figure  5,  with  the  white  area  denoting  the  admissible  region. 


EXPERIMENTAL  VALIDATION 

The  experimental  data  used  to  validate  this  work  is  provided  by  Professor  Schalk  Els  of  the  University  of  Pretoria, 
South  Africa.  The  tire  used  is  a  Continental  Conti-Trac  AT  238/85  R16.  Two  quasi-static  tests  (a  flat  plate  test  and  a 
cleat  test)  are  performed  that  produce  resultant  spindle  force,  F,  with  respect  to  the  tire  deflection,  e.  The  side  length 
of  the  square  tube  is  19  mm.  The  stiffness,  kg,  and  shape  parameters,  ary  and  a2,  were  tuned  to  approximate 
experimental  results  shown  in  Figure  5  and  Figure  6.  It  must  be  noted  that  a  more  formal  process  to  identify  the  best 
parameterization  of  the  model  is  under  development.  In  both  figures  the  experimental  results  are  shown  as  a  solid 
line  while  the  simulation  results  are  shown  as  a  dashed  line.  It  is  clear  from  Figure  5  that  the  nonlinear  stiffening 
characteristic  of  the  tire  is  captured  by  the  constraint  mode  tire  model  as  indicated  by  the  increased  slope  with 
increasing  deflection.  The  difference  between  the  simulated  and  experimental  forces  is  less  than  4%. 
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Figure  6.  Flat-plate  experimental  and  simulated  force,  F,  with  respect  to  the  deflection,  e 

The  cleat  test  results  shown  in  Figure  6  show  reasonable  agreement  throughout  the  range  of  deflection.  It  is  clear 
that  a  discontinuity  in  stiffness  occurs  when  the  deflection  exceeds  approximately  5cm,  dividing  the  results  into  two 
regions.  In  the  first  region  (displacements  up  to  5cm)  the  tire  is  suspended  by  the  cleat  and  the  contact  area  related 
to  force  generation  is  constrained  to  the  cleat  surface  only;  the  resulting  stiffness  is  relatively  small.  This  is  shown 
conceptually  in  the  upper  inset  picture  in  Figure  6.  Here  the  difference  between  the  experimental  and  simulated 
force  data  is  less  than  7%.  In  the  second  region  the  tire  tread  touches  the  surface  below  the  cleat  and  the  total 
contact  area  increases,  as  shown  in  the  lower  inset  picture.  This  discontinuity  in  stiffness  is  captured  by  the 
constraint  mode  tire  model.  Note  that  although  a  simple  linear  model  is  used  to  predict  the  radial  deflection  of  the 
tire  surface,  both  of  the  nonlinear  effects  shown  in  Figure  5  and  Figure  6  are  captured. 
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Figure  6.  Cleat  experimental  and  simulated  force,  F,  with  respect  to  the  deflection,  e 


DISCUSSION 

The  two  main  contributions  of  this  work  are  the  description  of  the  linear  planar  tire  model  in  terms  of  a  single 
stiffness  parameter  and  two  shape  parameters  and  the  development  of  an  admissible  region  for  the  shape  parameters. 
It  is  clear  that  the  particular  choice  of  tire  parameters  used  in  this  work  produces  reasonable  force  predictions  and 
accounts  for  nonlinear  changes  in  stiffness.  A  parameterization  process  that  more  exactly  replicates  the  experimental 
results  for  a  particular  tire  is  outside  the  scope  of  this  work,  but  is  the  focus  of  future  studies. 

It  is  envisioned  that  the  constraint  mode  tire  model  developed  in  this  work  can  serve  as  a  terrain  morphological 
filter.  That  is,  this  constraint  mode  tire  model  can  be  used  to  pre-filter  the  terrain  surface  once,  providing  the 
required  bridging  and  enveloping  properties  of  the  tire.  A  simpler  tire  model  (perhaps  a  simple  linear  point- 
follower)  could  be  used  in  the  iterative  design  process  to  provide  fast  yet  reliable  spindle  force  prediction  for  vehicle 
dynamic  simulation  and  reliability  evaluation.  Specifically,  an  effective  road  profile  could  be  estimated  given  an 
actual  terrain  profile  by  simulating  the  constraint  mode  tire  model  traveling  over  an  actual  terrain  profile  given  a 
constant  spindle  load.  It  is  hoped  that  a  small  sacrifice  in  accuracy  will  results  in  an  order  of  magnitude  increase  in 
computational  speed.  Ultimately,  the  tire  constraint  model  will  be  further  improved  so  that  it  can  be  applied  to  both 
deterministic  terrain  profiles  and  stochastic  terrain  profiles. 

CONCLUSIONS 

The  first  contribution  of  this  work  is  a  computationally  efficient,  planar,  quasi-static,  constraint  mode  tire  model  that 
captures  the  deformed  tire  shape.  The  required  model  parameters  are  reduced  to  two  non-dimensional  shape 
parameters  and  an  overall  stiffness  factor.  The  second  contribution  is  an  admissible  domain  for  the  shape  parameters 
that  are  sufficient  to  produce  a  realistic  tire  deflection.  Although  the  focus  of  the  work  is  to  develop  a  simple  linear 
model  to  predict  the  tire  shape,  simulated  spindle  loads  are  compared  to  experimental  loads  for  a  flat  plate  test  and  a 
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cleat  test.  The  results  of  the  experimental  comparison  demonstrate  that  nonlinear  stiffness  effects  are  adequately 
captured  by  this  simple  linear  tire  model.  A  computationally  efficient,  planar  tire  model  that  accurately  predicts  the 
lower-frequency,  but  not  necessarily  low  amplitude,  tire  shape  has  been  developed.  This  model  strikes  a  balance 
between  heuristic  tire  models  (such  as  a  linear  point-follower)  that  lack  the  fidelity  to  make  accurate  chassis  load 
predictions  and  computationally  intensive  models  that  cannot  provide  timely  predictions.  It  is  hoped  that  this  tire 
model  is  used  as  an  integral  part  of  the  computationally  efficient  and  accurate  vehicle  dynamic  simulations  that  are 
critical  throughout  the  iterative  design  process. 
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